Low Temperature Transport through a Quantum Dot: 
The Anderson Model out of Equihbrium 
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abstract 

^ The infinite-t/ Anderson model is applied to non-equilibrium transport through 
'"^ quantum dot containing two spin levels weakly coupled to two leads. At low 
SeJmperat ures, the Kondo peak in the equilibrium density of states is split upon the 
. application of a voltage bias. The split peaks, one at the chemical potential of each 
^ead, are suppressed by non-equilibrium dissipation. In a magnetic field, the Kondo 
peaks shift away from the chemical potentials by the Zeeman energy, leading to 
an observable peak in the differential conductance when the non-equilibrium bias 
equals the Zeeman energy. 

PACS numbers: 72.15. Qm 73.40.Gk 73.20.Dx 73.50.Fq 
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The behavior of an atomic impurity coupled to conduction electrons has become 
one of the paradigms of condensed matter physics. Competition between on-site 
Coulomb interaction and band hybridization produces the Kondo effect: a crossover 
from weak to strong coupling between the localized and band electrons below the 
Kondo temperature, Tj^ . The study of the Kondo effect has been limited, however, 
by the nature of the impurity system. Since it is a daunting task to drive the host 
metal out of equilibrium, it is the equilibrium properties of Kondo impurities that 
have been explored.^ 

In this paper we address a new Kondo system in which non- equilibrium is rou- 
tinely achieved, namely a semiconductor quantum dot weakly coupled to its leads. 
It is already evident that Anderson's modeP for a Kondo impurity - discrete, in- 
teracting levels coupled to a band - also describes quantum dots. Experimen- 
tally, the discrete spectrum of a single dot has been probed by transport^"^ and 
capacitance^ spectroscopy, while the strong on-site Coulomb interaction is observed 
in Coulomb-blockade conductance oscillations. ^'^'^ Theoretically, Anderson's model 
has provided an excellent description of these experiments both in equilibrium,^'^ 
and non-equilibrium.^'^ However, it is only the high temperature regime that has 
been explored experimentally, while it is below that the Kondo effect emerges. 

Since the Anderson Hamiltonian describes the quantum dot, at low temperatures 
the dot must behave as a Kondo impurity. In fact, Glazman & Raikh^^ and Ng & 
Lee^^ have argued that at zero-temperature equilibrium the Kondo resonance in the 
density of states of spin-degenerate levels will always lead to perfect transparency 
of the quantum dot at the Fermi energy. In contrast, above the Kondo tempera- 
ture, resonant tunneling occurs only at a discrete set of Fermi energies. Thus the 
Kondo effect will have a striking experimental signature in low-temperature trans- 
port through a quantum dot. Furthermore, in the quantum dot system the leads 
coupled to the dot are easily biased to non-equilibrium and the dot potential can 
be swept continuously with a gate. Thus new physical questions which were not 
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relevant to magnetic impurities can be raised. In particular, what happens to the 
Kondo effect out of equilibrium Since transport measurements on single quantum 
dots require significant applied bias, this question is of immediate importance. 

In this Letter we combine several approaches (non-crossing approximation,^^ 
equations of motion, perturbation theory, variational wavefunction calculation^^) 
to present a consistent picture of low-temperature, non-equilibrium transport 
through a quantum dot. For spin-degenerate levels at equilibrium, the Kondo peak^^ 
in the density of states at the chemical potential (Fig. la) leads to resonant trans- 
mission through the dot.^^'^^ A voltage bias between the left and right leads causes 
the Kondo peak to split, leaving a peak in the density of states at the chemical po- 
tential of each lead (Fig. lb). We find that the amphtudes of the split Kondo peaks 
are suppressed by a finite non-equilibrium lifetime. This lifetime results from dissi- 
pative transitions in which electrons are transferred from the high chemical potential 
lead to the low chemical potential one. As the voltage bias is increased, the hfetime 
decreases, resulting in a suppression of the Kondo peaks and thus a suppression of 
the differential conductance (Fig. 2a). Upon application of a magnetic field, the 
Kondo peaks shift away from the chemical potentials by the Zeeman splitting, but 
in opposite directions for each spin (Figs. Ic and Id). Interestingly, therefore, when 
the chemical potential splitting equals the Zeeman splitting, a Kondo peak shifted 
away from one chemical potential crosses the other chemical potential. We predict 
an observable peak in the differential conductance at this crossing^^ (Fig- 2b). 

We model the quantum dot and its leads by the Anderson Hamiltonian^ 

H= ekaCt^Cka + J2e^c+c^ + Un^n^+ J2 {VkaCtaCa + h.c), (1) 

where cl^{ck(j) creates (destroys) an electron with momentum k and spin cr in one of 
the two leads, and cj(ccr) creates (destroys) a spin-u electron on the quantum dot. 
Since we are interested in temperatures smaller than the orbital level spacing in the 
quantum dot, we consider only a single pair of levels on the dot with energies = 
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eo+Ae/2 and = eo— Ae/2. The third term in (|l]) describes the Coulomb interaction 
between the two locahzed spins which we take to forbid double occupancy^^ {U —>■ 
oo), while the fourth term describes the hopping between the leads and the dot. 

Our aim is to calculate the current through the dot, J, which for the case 
of proportionate couplings to the leads, T^^u) = \T^{uj), where r^'^'(ci;) = 
'^'^Y^keuR) \^ka\'^S{^ — e/ccr), can be expressed^° in terms of the density of states. 



In Eq. (I), T„{u) = T^^{uj)T^{uj) / [r^(cj) + T^{uj)], and G^{uj) is the Fourier trans- 
form of the retarded Green function, G^(t) = —iQ(t){{c„(t),c'^{0)}). 



In order to calculate the Green function G^{uj) we use both the non-crossing 
approximation^^ and an equations-of-motion method. The non-crossing approx- 
imation is based on an exact mapping of the infinite- [/ Anderson Hamiltonian (|1|) 
onto a slave-boson Hamiltonian. If vertex corrections are neglected, the propagators 
for the boson and the fermion degrees of freedom, which correspond, respectively, to 
the propagators for the empty site and a singly occupied site, obey a set of coupled 
integral equations. Numerical solution of these equations has been very useful in ob- 
taining quantitative results for the equilibrium system, including the occupations 
of the two spin-states in the presence of a magnetic field. In this work we have 
generalized the non-crossing approximation to non-equilibrium to produce densities 
of states, occupations, and the nonlinear current (Q). However, as a large spin- 
degeneracy (large N) technique, the non-crossing approximation produces a Kondo 
peak even for the non-interacting system (A^ = 1). Consequently, for = 2 in a 
magnetic field, it give rise to spurious peaks in the density of states at the chemical 
potentials. Therefore, an equations-of-motion method was employed to comple- 
ment the non-crossing approximation and isolate its shortcomings. This method 
corresponds to a resummation of low-order hopping processes and cannot produce a 
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quantitative description of the Kondo effect. Nevertheless, this method is known^^ 
to give the right quahtative behavior at low temperatures. More importantly in the 
present context, since the equations-of-motion method is exact for A'" = 1, it gives 
rise only to the proper Kondo peaks (as identified by perturbation theory^^). 

The equations-of-motion method consists of differentiating the Green function 
G^{t) with respect to time, thereby generating higher-order Green functions which 
eventually have to be approximated to close the equation for G^{t). The procedure 
we employ here is the same as the one used in Ref. 8, which in the infinite-C/ limit 
leads to 



with 



and 



^oaH = E '^'^'1 - , (4) 



where fL/ni^) is the Fermi distribution in the left/right lead and is the 
intermediate-state lifetime. G^{u!) has an overall amplitude proportional to 1 — (n^) , 
where (ria-) is the occupation of the other spin-state. Quantitative calculation of the 
occupations is beyond the scope of the equations of motion in the present approxi- 
mation scheme. Accordingly, we use the occupations resulting from the non-crossing 
approximation, which are known to be quantitatively reliable in equilibrium.^^ 

Within the equations-of-motion scheme, the Kondo peak for spin a results from 
the self-energy, Eio.(a;), due to virtual intermediate states in which the site is occu- 
pied by an electron of opposite spin, a. The remaining self-energy, Eo(t(<^), is the 
exact self-energy for the non-interacting case. Because of the sharp Fermi surfaces 
at low temperature, Re{Eio.(a;)} grows logarithmically at a; = //^^^i Ae, giving rise 



to peaks in the density of states, —^lmG^{u;), near those energies. The peaks for 
the high- lying spin (low-lying spin) appear near cu = fir^ n + Ae {uj = fii^ n — Ae). 
At zero-field and zero-temperature equilibrium, the intermediate states giving rise 
to Eio.(a;) have an infinite lifetime, and the true peak in the density of states has 
an amplitude corresponding to the unitarity limit. Once either a voltage bias or a 
magnetic field is applied these intermediate states acquire a finite lifetime, r^, which 
cuts off the logarithmic divergence of Re {Eio.(a;)}, resulting in a suppression of the 
peak amplitudes. The lifetime, Ta, of spin a can be calculated using second-order 
perturbation theory, yielding 

1 1 1 roo 

Ta n A=L,R A,B=L,R •'-oo 



1 1 

+ 



Sa-e + irj)'^ {ea-e-ir}y_ 
r^(e)rf,(e -ea + e^O (1 " /^(e)) /^(e - + e,0" • 



For a deep level at zero temperature and for constant F this simplifies to 

— ^TTT E ^a^a'^il^B- I^A + ea-ea')-, r, 7) 

Ta 27rft [flA- ea){HB - €a') 

a' 

which explicitly shows that the lifetime is non-zero only for finite bias or finite 
magnetic field. 

In Fig. 1, we plot the density of states for two spins symmetrically coupled to 
two leads, consisting of Lorentzian bands of width 2W, so that F^(a;) = F^(a;) = 
rW'^/2(u;'^ + W'^), with F = 1 and W — 100. Results are shown for the non-crossing 
approximation (dashed lines), which is reliable for zero magnetic field, and for the 
equations-of- motion method (continuous lines), which has the correct Kondo peak 
energies for all magnetic fields. In equilibrium and zero magnetic field, the density 
of states exhibits a single peak at the Fermi level as expected^^ (Fig- la). As the 
chemical potentials split, the Kondo peak also splits, giving rise to a suppressed 
Kondo peak at each chemical potential (Fig. lb). Upon the application of a mag- 
netic field, the densities of states for the two spins become different and the Kondo 
peaks shift away from the chemical potentials by the Zeeman splitting (Ae = 0.2 in 
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Figs. Ic and Id). The peaks move up in energy for the high-lying spin (Fig. Ic) 
and down in energy for the low-lying one (Fig. Id). 

The main conclusion of Fig. 1 is the emergence of new energy scales, not present 
in equilibrium. The Kondo peak in the equilibrium density of states splits out 
of equilibrium to two peaks spaced by the chemical potential difference A/i, and 
suppressed from equilibrium by the finite dissipative lifetime To-. In Fig. 1(b), the 
lifetime broadening, ^/to-, is about the same as the temperature. It is apparent, 
however, from Figs. 1 (c) and (d), that neither the non-crossing approximation nor 
the equations of motion quantitatively determine the Kondo peaks at finite magnetic 
field. For this case we use the equations-of-motion result since it provides a good 
estimate of the Kondo peak positions (by comparison with perturbation theory). 
To understand the shift of the Kondo peaks with magnetic field, it is helpful to 
recall how the peaks in the density of states derive from the eigenstates of the 
system. At T = 0, G'^{t) involves transitions from the A^-particle ground state to 
all possible + 1 or — 1 states. At B = the correlated ground state has a finite 
amplitude to have an empty site, and thus c+ (c^) can generate transitions from the 
A^-particle ground state to the ground state with one more (one less) electron. Since, 
by definition, the ground state energies differ by the chemical potential, the density 
of states includes a Kondo peak at the chemical potential. Within a variational 
calculation,^^ we find that at finite magnetic field the ground state is polarized, and 
adding or removing an electron produces no overlap with the new ground state. 
However, there is a correlated excited state of opposite polarization which can be 
reached, and which consequently gives rise to a peak in the density of states, shifted 
by the difference in energy between polarization states, i.e. the Zeeman energy. 

The current follows immediately from the densities of states (Q). In particular, 
the zero-temperature current is the integrated density of states between the two 
chemical potentials, weighted by the coupling to the leads Ta-i^uj). At zero mag- 
netic field, therefore, the Kondo peak at the Fermi energy gives rise to a linear- 
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response conductance of 2e^/h for symmetric barriers, corresponding to perfect 
resonant transmission through the quantum dot.^^'^^ As the bias is increased the 
differential conductance falls rapidly (Fig. 2a). This occurs firstly because the 
differential conductance due to a peak in the density of states must fall off once A// 
exceeds the peak width, and secondly because the decreasing dissipative lifetime 
suppresses the peak amplitudes. Since the peaks in the density of states persist 
until the temperature is roughly one-tenth the coupling to the leads, F, the peak 
in the differential conductance is observable well above the Kondo temperature, 
(Fig. 2a, continuous line). 

In a finite magnetic field the Kondo peaks are shifted away from the chemical 
potential so they contribute very little to the conductance in linear response. As 
the bias is increased, however, the current carrying region between the chemical 
potentials grows, until at A// = Ae, it reaches one Kondo peak in the density of 
states of each spin (see inset of Fig. 2b). In Fig. 2b, one therefore sees peaks in 
the differential conductance at A/x = Ae (continuous line). In fact, by comparison 
with the non-crossing approximation (Fig. 1), we expect the equations of motion 
to underestimate the full strength of these peaks. Experimentally, observation of 
peaks in the differential conductance at A/x = Ae would provide a "smoking gun" 
for the presence of Kondo physics in transport through a quantum dot. 

In this work, we addressed the non-equilibrium behavior of Anderson's model for 
a magnetic impurity. Experimentally, the model describes low-temperature trans- 
port through a quantum dot, where non-equilibrium is readily accessible. We have 
shown that new energy scales emerge in non-equilibrium. Specifically, the difference 
in chemical potentials A// and the inverse dissipation time h/r^ lead, respectively, 
to splitting and suppression of the Kondo resonances in the density of states. Our 
results have led to a novel experimental prediction — when the Zeeman splitting of 
the spins, Ae, equals the applied bias, A/i,, there will be a peak in the differential 
conductance, provided these energies are smaller than the coupling to the leads, F, 
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and smaller than the depth of the levels, //i,^ — e^- Importantly, this signature of 
the Kondo effect persists, for a wide range of parameters, to temperatures ~ r/20, 
which may be magnitudes larger than the Kondo temperature. We hope that this 
work will encourage further efforts, both experimental and theoretical, to probe the 
non-equilibrium physics of interacting quantum systems. 

Acknowledgments: We thank W. Kohn, T. K. Ng, and N. Read for valuable 
discussions. Work at UCSB was supported by NSF grant no. NSF-DMR90-01502 
and by the NSF Science and Technology Center for Quantized Electronic Struc- 
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Figure Captions: 

(1) Density of states for an Anderson impurity symmetrically coupled to two leads with 
chemical potentials /jLl and /xr(= 0) and Lorentzian bandwidth 2W, from the equations-of- 
motion method (continuous line) and the non-crossing approximation (dashed line). The 
impurity has two spin states with energies e-f and and an on-site interaction U ^ oo. 
All energies are in units of the total coupling to the leads, T. The band width is W = 100 
and the temperature is T = 0.005, roughly a factor of two lower than the magnetization 
Kondo temperature [Ref. 14]. (a) The equilibrium {/il = 0) density of states at zero 
magnetic field = = —2.0. The density of states exhibits a single peak at the Fermi 
level [Ref. 17]. (b) The non-equilibrium (/x^ = 0.3) density of states at zero magnetic field 

= El = —2.0. There is a suppressed Kondo peak at each chemical potential. (c),(d) 
The non-equilibrium (/xl = 0.3) density of states for spin up (c) and spin down (d) at 
finite magnetic field = —1.9, = —2.1. The Kondo peaks shift away from the chemical 
potentials by the Zeeman splitting Ae = 0.2; the shift is up in energy for the up spin and 
down in energy for the down spin. 

(2) Differential conductance, edJ/dAn, with = 0, vs. applied bias, (a) Zero magnetic 

field differential conductance via the non-crossing approximation, (b) Differential conduc- 
tance at the finite magnetic field, Ae = 0.2, used in Figs. 1 (c) and (d), via equations of 
motion. As shown in the inset, when the chemical potential difference, A/x, reaches the 
Zeeman splitting, Ae, the Kondo peaks in the density of states enter the region between 
the chemical potentials, giving rise to a peak in the differential conductance. 
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